function argmax = s_dot_max(dV_s_over_ds)
global delta
global gamma_x
global gamma_s
global alpha
global c_s_1
global c_s_2
global x
global s
global s_dot_bar

state_dependent_term = max(gamma_s-s,0)-max(gamma_x-x,0);

argmax = min(max(0,((dV_s_over_ds-state_dependent_term-c_s_1)/(c_s_2*alpha))).^(1/(alpha-1))-delta,s_dot_bar);

end

